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In January 1964 an investigation was begun to investigate the 
validity of classical rigid aircraft analytical techniques on predicting 
the stick fixed and dynamic longitudinal stability of the Goodyear YAO-3G 
Inflatoplane. In addition, a study was made on the effects of fuselage 
flexibility on the dynamic response in the phugoid mode and the movement 
of neutral and maneuver points. 

It was found that rigid techniques were valid for dynamic response 
predictions and somewhat less valid for static predictions. However, if 
the airframe pressure was maintained at the designed level static pre- 
dictions were in no way adverse. If the airframe is available for static 
bending test, correction factors may be obtained that when applied to 
rigid values of Cmj, result in more accurate predictions of both dynamic 
and static stability. 

Fuselage rigidity was varied by changing internal airframe pres- 
sure from 5 psig to 7 psig. The effect of fuselage flexibility was to re- 
duce damping and increase the frequency of the phugoid mode toa slight 
extent. The neutral and maneuver points both moved forward; however, 


it was indicated that iets moved further forward than No 








INTRODUCTION 


The Goodyear YAO-3G is an aircraft that was developed in the 
1950's and is no longer in production. The aircraft is unique in that it 
is of inflatable rubber-fabric construction. 

In October of 1963 the Department of Aerospace and Mechanical 
Sciences of Princeton University obtained one of these aircraft on loan 
from the U. S. Army. It was determined that the unique features of the 
aircraft could be put to good uSe in an investigation of first order aero- 
elastic effects due to structural bending. The investigation would be 
simplified by the low aerodynamic pressure required to manifest such 
bending. In addition, the rigidity of the aircraft could be varied, with- 
in limits, by controlling the internal pressure of the airframe. 

With these possibilities in mind this investigation was started in 
January of 1964. The first flight test was conducted on January 3, 1964 
and the tests continued through early May of 1964. 

The purpose of the investigation was twofold. First to investigate 
the validity of classical, rigid aircraft, analytical techniques in pre- 
dicting static and dynamic longitudinal stability, and secondly to investi- 
gate the effects of fuselage and wing bending on the stick fixed maneuver 
and neutral points and the dynamic modes of longitudinal motion of this 
aircraft. 

The investigation consisted of analytical developments consider- 
ing the aircraft first rigid and then semi-rigid and comparing the 


analytical analysis with flight test results. 








EQUIPMENT 


The Aircraft 


The YAO-3G is a single-placed high-wing monoplane of inflatable 
fabric construction (Figures 1 through 6 and Table I). The landing gear 
consists of a single fixed wheel with wing tip and tail skids. The aircraft 
is powered by a single engine driving a fixed pitch propeller. The engine 
is a Nelson Model H63-A which develops a maximum of 42 brake horse- 
power at 4000 RPM at sea level. The engine operates on the two-cycle 
principle, has a single magneto ignition, and air pump for controlling 
airframe pressure. 

Other than being inflatable, the airframe is of conventional form 
consisting of a straight constant chord wing (NACA 0015), a fuselage with 
partially enclosed cockpit, and empennage with vertical and horizontal 
stabilizer and control surfaces. A conventional stick and rudder pedal 
control system is used for flight and manual control is provided for the 
engine mounted air pump and airframe pressure relief valve. 

The rigidity of the aircraft is maintained by hoop stresses result- 
ing from internal pressures varying from 4.8 to 7.0 psig. The shape of 
the non-cylindrical parts is maintained by drop stich threads running 
from one surface to the other as shown in Figure 6. The internal pres- 
sure changes with altitude and is maintained at the desired level by 
manipulation of the pump control and airframe pressure relief valve. 

The instrumentation is quite basic consisting of an air speed 
indicator, altimeter, accelerometer, magnetic compass, air pressure 
gage, engine tachometer, and cylinder head temperature indicator 
(Figure 7). The electrical system consists of a dry cell battery powered 


warning light for the air pressure gage. 








The basic weight of the aircraft was 315 pounds including test 
instrumentation and parachute. Gross weight varied from 500 to 


600 pounds depending on pilot and fuel. 


F lig ht Test Instrumentation 


The instrumentation requirements were to provide the necessary 
data to perform static and dynamic longitudinal stability analyses. This 
required a method of measuring and recording elevator deflection, nor- 
mal acceleration, airspeed, angle of attack, and pitch rate. In addition, 
a movie camera would be desirable to study fuselage bending during the 
phugoid motion, and voice communications requirements dictated self- 
contained radio equipment. 

The instrumentation problems were somewhat different than would 
be expected of a more conventional aircraft. The space available for in- 
strumentation was a small area behind the pilot on the cockpit floor. This 
space had to contain the major instrumentation components and power sup- 
ply. Weight was also a problem in that a 200 pound pilot plus parachute 
and nominal fuel put the aircraft at near maximum gross weight. In 
addition, attachment points were limited and some method had to be de= 
vised to attach and remove instrumentation components without damage 
to the fabric airframe. 

A SFIM A20 Flight Recorder was chosen to record data. This 
instrument is 6.5 x 4.7 x 3.9 inches and weighs 4.6 pounds including air- 
speed and acceleration transducers. A matched rate gyro 3 x 3 x 3 inches 
weighing 3 pounds was used to record pitch rate. Transducers in the re- 
corder convert inputs to mirror deflections. The mirrors reflect light, 
from a single source, to light sensitive paper. Five input channels may 
be recorded on the 60 mm paper. The instrument and associated equip- 


ment draws 2 amps at 28 volts D.C. 






The A2O can be set to run at tape speeds of 1 mm or 5 mm per 
second and holds approximately 25 feet of recording paper. To compen- 
sate for variations in tape speed with supply voltage, an internal clock 
mechanism is provided which projects time marks on the trace. In addi- 
tion, a trace marking device is incorporated. 

The recorder transducers consisted of an airspeed indicator that 
required static and dynamic pressure inputs, a self-contained accelero- 
meter, two ratiometers, anda galvanometer for potentiometer inputs. 
The ratiometers compensate for variations in supply voltage and were 
used for pitch rate and elevator deflection inputs. The galvanometer is 
not insensitive to supply voltage and was used for the angle of attack in- 
put in that a qualitative value of angle of attack would be sufficient for 
the dynamic test. A more complete description of the SFIM A20 may be 
found in Reference 1. 

The power supply consisted of 24 nickel cadmium batteries witha 
nominal output of 29 volts and capacity of 4 amp hours. The power sup- 
ply was made up in two packs of 12 batteries each in order to better 
utilize space requirements. 

The mounting problem was solved by the use of Permacel P69 
fabric tape. This tape was found invaluable for manufacturing attach- 
ment points. It was found to adhere well through the various tempera- 
ture and dynamic pressure conditions encounted during flight test. 

Figure 8 shows the installation of the A20 recorder, rate gyro, 
and power supply. 

In order to accurately measure elevator deflection, it was neces - 
sary to put a potentiometer on the elevator itself. Due to fuselage bend- 
ing, the elevator could deflect with the stick fixed and cable measurement 


was not considered to be of sufficient accuracy. The potentiometer 






mounting is shown in Figure 9. The potentiometer leads were taped to 
the fuselage right side. Due to the low dynamic pressure of the aircraft, 
external mountings were not a problem. 

The mounting of the angle of attack vane presented a problem as 
to what part of the aircraft would be rigid enough to use as a reference. 
A typical wing mounted boom was out of the question due to wing twist 
and flexibility. As a compromise solution the boom was mounted by a 
two-foot bracket attached to the landing gear frame. The gear frame is 
the most rigid part of the aircraft; however, in order to preserve the 
rigidity of the mounting bracket, the boom had to be shorter than neces- 
sary to totally escape from wing interference. However, it was con- 
sidered that the amount of wing interference encountered would not ad- 
versely affect angle of attack dynamic response. Although primarily 
designed as a mounting for the angle of attack vane, the mounting bracket 
seemed to be the ideal location for the movie camera (a surplus gun 
camera), also. Later on it was also used to mount a venturi tube. This 
venturi was necessary to amplify static pressure for the recorder air- 
speed transducer. The transducer was found to be insensitive at the 
low airspeeds encountered when attached to a simple pilot static tube. 
Figure 10 shows the angle of attack, camera, and venturi mounting. 

A control panel for the recorder was mounted on the instrument 
panel in place of the magnetic compass (Figure 7). This panel contained 
a master switch to control power to the recorder motor and lamp, a 
speed switch that controlled the recorder motor, an event marker switch, 
and a five position selector switch and meter. The selector switch 
allowed the meter to read either elevator position, angle of attack, pitch 
rate, battery voltage, or off. This selector arrangement proved in- 
valuable in that recorder inputs could be checked in flight. In addition, 


a switch panel was mounted on the right cockpit wall. This panel 






contained angle of attack, pitch rate, and camera switches. The position 

of the angle of attack and pitch rate switches determined whether these 
values were being recorded. As these values were not necessary for the 
static test, their presence made trace reading more difficult. The switch 
panel allows static and dynamic tests to be conducted on the same flight 
without excessive traces on the tape. When the camera switch was on, 

the recorder motor switch also controlled the camera. An event marker 
switch was also mounted on the top of the pilot's control stick. Figure // 
is a wiring diagram of the complete instrumentation package and Figure 12 
is a schematic of the recorder system. 

The voice communication system consisted of separate self-con- 
tained VHF transmitter and receiver. The transmitter was a Navguide 
Corporation Model AT -103 which was 5 x 7 x 2 inches and weighed 
5 pounds including batteries. The Navguide is a crystal tuned 5 tube 
VHF transmitter which operates on dry cells. The filaments do not 
light until the microphone button is depressed. It was found that the 
filament battery would not sustain semi-continuous operation necessary 
during flight test. To aleviate this the recorder power supply was put 
in parallel with the filament battery after being dropped down to 14 volts 
with a 4000 ohm, 5 watt resistor. The transmitter was mounted behind 
the pilot just under the leading edge of the wing. The transmitter antenna 
was taped to the left outer cockpit wall (Figure 3). The receiver was a 
Montery Electronics "Partner" continuous tuning VHF receiver 2x4x5 
inches weighing 13 pounds including mercury batteries. The receiver 
was mounted on the left inner cockpit wall (Figure 7). The receiver 
antenna was taped to the left outer cockpit wall (Figure 3). 

As the aircraft has no metal skin there was no ground plane effect 


for the antenna mountings and the receiver would not operate under these 






conditions. A ground plane was constructed from 1/16 inch welding rod 
and taped to the underside of the cockpit. This arrangement proved quite 
satisfactory for two way communication up to 5 miles. As this range 
was satisfactory for flight test no ground plane was used for the trans- 
mitter antenna. 

The total weight of the instrumentation including all mountings 
amounted to 33.5 pounds. The SFIM recorder operated very well 
throughout the flight test. Its light weight, compact size, and simple 


power supply requirement solved the basic instrumentation problem. 






ANALYTICAL PROCEDURES 


Longitudinal Stick-fixed Stability 


The longitudinal stability is determined by considering the motion 
of the aircraft in horizontal and vertical translation and in pitch. Symmetry 
of the aircraft precludes any coupling between longitudinal and lateral 
motions. The wind-axis is assumed to be parallel to the principal longitu- 
dinal axis in the vertical plane of symmetry, thus eliminating the products 
of inertia. The resulting perturbation equations of motion are, from 


Reference 2, 
Vv + = TT" V+ a= cs)) + pA@ = 
A (Dy, TY)A (D  g)he gA 0 


L ne . 


V Qa 

= ween - A@ = 

v AV + Aa y Oe A 0 
° © 


-M)AV -M ha ~-M da + A@ - Mg 48 = Ava 


In this evaluation the aircraft will be considered as a rigid body 
except where full-scale wind tunnel test data is available. Using the air- 
craft specifications of Table I and references as noted, the stability deriva- 
tives may be evaluated as follows. The flight conditions used for evalua- 
tion will be: 

W = 515.2 Ibs. 
c.g.= 29% mac 
h = sea level 


V = 50 knots = 84.45 fps 
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The lift coefficient is determined from: 


ee Zils (Z)(50522) 

= Rela 7 RASS Dna OE nana 
L pv’sS (.00238)(84,. 45)* (124. 3) 
C,, = + 488 


The drag coefficient is determined from: 


2 
G = G + em 
Due Do eTAR 


With Cn, = .1l from Reference 3, and a reasonable value of the aircraft 


efficiency factor, e, takenas .75, the drag coefficient becomes: 


5 (. 488)? 
Cp 7) eG 1414044) 


Gey = 23 


The aircraft non-dimensional mass based on wing chord is: 


m 16 
[e-  CoeiensiGen 22 


The effect on the forces and moments due to a change in forward 


speed with angle of attack and throttle fixed are stated below. Drag 


changes with speed by: 








bi 
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Dre, 1915 


The thrust versus velocity curve shown in Reference 3, gives a change of 


thrust with velocity of -. 711 lbs/ft/ sec. 


1 dT -, 711 
Voemdv 16 


Tp =, 04445 
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The change of lift with a change of velocity is expressible by: 


nae 2g _ (2)(32.2) 
Ve ~ Vv" ~ (84, 45)" 
fe 

We 
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With the engine lever arm 3.14' above the center of gravity and assuming 
the thrust line parallel to the fuselage reference line, the change in pitch- 


ing moment with speed is: 


_ 1 8M_1 aT, _ (-.711)(-3.14) 
Vo lyeoy Iveov eee 


= ,0080 
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The changes that occur in the forces and moments with a change in 
angle of attack are normally an increase in lift, anrincrease in drag and a 
negative pitching moment. The lift curve slope of the aircraft at 50 knots 


as determined from full-scale wind tunnel tests, Reference 4, is: 


oe = 4.18 /rad 
a 


The drag change with angle of attack is represented by: 


p =28_ ¢ = (232.2)(4.18) 
a etAR Ly (. 75)(3. 14)(4. 44) 


D = 26,0 
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A change in angle of attack effects the lift by: 


x: _ (32.2)(4.18) 
V.. C_V._. (.488)(84. 45) 
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In determining the angle of attack effect on the pitching moment the con- 
tributions of the fuselage and tail must be considered first. The fuselage 


contribution may be estimated from Reference 5 by: 
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Where L. is the overall fuselage length, Me is the maximum width of the 


fuselage, and K, is an emperical factor determined from experimental 


if 
evidence. 


(. 012)(2. 25)°(19. 6) 
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Cc 
An angle of attack change at the tail produces a moment change as follows: 


ma . 
The tail efficiency, ale = = is felt to be about 1.3 due to slipstream effect. 
This value was verified experimentally as described in the test procedures 
section and found to be a reasonable assumption. The lift curve slope of 


the tail can be estimated from Reference 6 to be, Ps 257/ rade 


Gang, = ~(2- TM. 3). 399) = -1. 397 


A lag in the downwash reaching the tail may be found from the emperically 


determined formula: (Reference 14) 
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Where y is the vertical distance of the tail quarter chord from the wing 


chord line (+ above). 
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Cma = (4.18)(. 04) + .0184 - 1.397 (1 - . 42) 
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M_ = -12.61 
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There is also a moment change due to rate of change of angle of attack 


given by: 
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The principal effect of pitch is manifested in a damping of the 


moment as follows: 
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The elevator control power may be expressed by: 
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where Le the control effectiveness parameter, is .6 from Figure 5-33 of 


Reference 7. 


_ (32.2)(. 6)(-1. 397) 
§  (.488)(5.31)(. 615) 


-16. 82 


< 
th 








16. 


Using the stability derivatives above, the following determinant of 


coefficients is obtained: 


AV Aa A@ 
S + 230 -6.2 32.2 
. 00903 53 265 -s 0 
-. 00806 1.7465 +12.61 s(s+4.16) 


From this the characteristic equation for the rigid aircraft is: 
A* +9. 407A + 28.4124A% + 7.1315A + 4.514 = 0 


Solution of the equation yields phugoid and short period characteristics of: 


Phugoid Short Period 
C= .245 € = 3897 
on = 414 on=5.128 
TS = 7.05 sec Tte=) . 15igsec 
Period = 15.7 sec Period = 1.22 sec 
we = .401 Wc = 2,28 


Static Stability 


The neutral point and maneuver point for the rigid aircraft can be 
determined from the stability derivatives once they have been converted 


to non-dimensional form. 
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The resulting stick-fixed neutral and maneuver points are: 


Ni = 47.5% mac 


N 
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59.9% mac 


The slope of the neutral point curve may be determined by differentiating 
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With the elevator effectiveness estimated as, i = .6, the slope of the 
neutral point curve is: 
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The slope of the maneuver point curve may likewise be deter- 


mined by differentiating —-— with respect to a 
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The slope of the maneuver point curve for the rigid case is: 
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Analog Computer and Root Locus 


Due to the flexibility of this aircraft, the vertical bending of the 
rear portion of the fuselage will manifest itself in a change in the hori- 
zontal stabilizer incidence, i, This change can be considered as a change 
in the derivative ome This effect has serious aerodynamic consequences 
because all the pitching moment stability derivatives, with the exception 


of speed stability, are a function of C,,,. - 
i 
t 








20, 


The effect of wing bending predominately manifests itself in a change 
in ale . In this investigation a varying oie will be used. 
In order to observe the result of varying tail incidence, ee was 


factored out of the moment equation. 
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This leaves the moment equation in a suitable form for both root locus and 


analog simulation. 
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The root locus equation for varying Cmi, was derived from the 


following determinant where . = a 
mu 
+ 4 - 
Mi Vv 
° oO 
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Ly s + vy aS 
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Upon substituting into the above determinant for all values except Cy and 


Cm; , the following determinant was formed. _ 
t 
sat 236 (.0736 C, + .3815) + 2516 N5 
@ 
aod st - 7oIeC ys, - s 
Q 
- .681 {.814C, - . 3625 s(s -2.977 Cy, ) 
a t 


= Cm, [1.251 6 + 14,693) 
A 
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This produced the following root locus equation for a varying C,,. 





1t 
4.238|Cmn, [s° +(2. 9943+. 5486Cy, )s* +(. 9406+. 0901 C, )s+. 8008 ] 
t Qa a 
ne ee eS ee Bl 
4 3 2 
+(.236+. 781 +(-. 0722 -.6777 +(=-. 0 -.14 +(-.1052 - 
s'+(.236+.781Cy,_ )s + CL, )st(-. 0855-. 14Cy, )B+(-. 105 
-0292C, ) 
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CL, was not substituted in because the root locus will be drawn for 
a range of ci 's. This procedure was used because ue has a large effect 
on the dynamics of the aircraft. Also the variable ci obtained from 
Reference 3 could not be depended on since the configurations do not match 


exactly. The root locus will be drawn for a range of possible Clay 's, 





Cc = 3.0, 4.0, and 5.0 per radian. 
Lg 
C, =3.0 
le 
4.238] Cm, [s°+ 4.6401s°+ 1.2110s + . 8008] 
L. 
t 


s*+ 2.579s°- 2.1053s8°- .5055s - .1925 
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= -l 
(s + 3.1941)(s - .8591)(s + .122 +j .2351) 
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4.2381C py, (s+ 5.1887s°+ 1.3012s + . 8008) 





c 





B43. 3659 2920 7836" = .64555 = 222116 


4.238] Cin, 
1 





(s + 4.9588)(s + .1149 + 3.3852) 
it 


(s+ 470166)(s -.8900)(s + 117 £7 .2197) 


4.238 (Crm, (s°+5.7373s* +1.3913s + .8008) 





t 


s*+4.141s* - 3.4607s° - .7855s -.2507 


4.238 (seo colle )\(s°+ , 131 2) 2357) 








Cm. 
1 

t 

= «] 


(s + 4.8268)(s - .9143)(s + .1143 +j.2092) 


The root locus for varying omi, for the three cn *s is shown in 
Figure 16. Figure 17 shows the stable phugoid portion of Figure 16 with 
an expanded scale. Lines of constant emu are shown plotted on these 
curves. Figure 17 also shows lines of constant damping ratio, ¢. 

The analog simulation was based on the normal lift and drag 


equations plus the modified moment equation. The equations used in the 


analog simulation are shown on the following page. 
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Figure 13 shows the computer set-up used. Table II is a tabulation of the 
pots versus the quantity they set. The computer set-up was made in such 
a way as to keep C15 and Cmi, separate. This could be accomplished 
easily for Cmi, by taking advantage of the previously factored equations. 
However, the oie could not be so easily separated since it enters the 
lift, drag, and moment equations. 

The computer set-up was checked by comparing the analytic solu- 
tion with the computer solution. This Sanborn tape is shown in Figure 14. 


As can be seen in Figure 17, the computer solution plots very close to the 


analytic solution. 
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Analog computer runs were made for ous 's of 3.0, 4.0 and 5.0 
per radian. ou was varied from 2.0 down to the Cmi, for critical 
damping. Typical Sanborn tapes are shown in Figure 15. 

Agreement between the root locus and the analog computer was 
checked by making analog runs at Cus "s of 3.0, 4.0, and 5.0 per radian 
and with Cy; =1.0. The W's were taken directly from the tapes while 
the ¢'s were found by taking amplitude ratios from the tapes and applying 
them to Figure 18. These values are shown plotted on Figure 17. 


All analog runs are fora oF of .15 radians introduced as a pulse 


of about 3 seconds duration. 
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FLIGHT TEST PROCEDURES 


The flight test phase of the investigation began in January 1964 
and continued through early May 1964 consisting of 33 hours of flight time. 
These tests were conducted at the air facility at the Forrestal Research 
Center of Princeton University. This phase consisted of preflight calibra - 
tions, in-flight calibrations, stick fixed longitudinal static stability flight 
test, stick fixed longitudinal dynamic stability flight test, flight test de- 
termination of pitching moment due to tail angle of attack (Cmi,)» and 


elevator effectiveness Ee 


Preflight Calibrations 


As soon as the test instrumentation was installed the aircraft was 
weighed with various fuel loads to determine fuel weight and moment con- 
tributions. The aircraft was then weighed with no fuel and each of the 
three pilots. The fuel weight and moment were then added to pilot weight 
and moment to arrive at center of gravity locations and gross weight 
versus fuel load (Figure 19). 

It was found that three pilots weighing 170 lbs., 190 lbs, and 
210 lbs., respectively, produced center of gravity excursions from 
24.10% MAC to 30.00% MAC. An attempt was made to extend this range 
of center of gravity locations aft by adding weights to the tail; however, 
this introduced nonlinearities in fuselage bending that could not be ac- 
counted for in the static stability test. It was decided that witha large 
selection of flight test points the static stability curves could be ac- 
curately extrapolated with the center of gravity excursion provided by 
different pilots. 

The elevator motion was calibrated by attaching a protractor to 


the elevator and recording deflections obtained by selecting control 
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stick positions. This calibration was found to be different for 7 psi and 

5 psi internal pressures. As elevator calibration was critical for static 
stability tests, the elevator calibration was done every time flight test 
data was obtained. A possible explanation for the change of elevator 
calibration is that the string attached to the potentiometer could change 
length due to aging and ambient weather conditions. The use of wire would 
eliminate the problem, but was ruled out for flight safety reasons. A 
string was chosen instead of a wire to preclude the possible jamming of 
the control system by the instrumentation. 

The only calibration required of the angle of attack vane was to 
assure that its nominal position fell in the center of the recorder trace 
as no quantitative results were required of it. 

The recorder accelerometer was simply recorded at +1 and -1 ¢ 
in order to check the flight test calibration for a 2g range. 

In predicting Cm; for the analytical development it was necessary 
to choose a value for tail efficiency (n,)- It was felt that Me for this air- 
craft would be greater than 1. 0 because there are no fuselage effects to 
dissipate the wake other than that obtained by the downwash effect off the 
wing lowering it into the fuselage. Considering no downwash effect the 


following development was made: 


From the 


momentum theory: = ee 





(1) Matec SW ab wy 


where v is velocity induced by disc at the disc 
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(2) ee = V + 2v (downstream) 


Therefore mass flow = PAV +v) 


and from (2): 





V - V 
ied ss 
_ 2 
Therefore: 
PAY 
mass flow =—>— (V__ + V) 
2 ss 
and 
=(PpA 
Thrust = (9 D V_)4v 
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and from the analytical development made at test conditions, i.e., 


V = 54 knots IAS 
Drag = 129.4 lbs = Thrust 


Therefore: 


(2)(129. 4) 


_ f bento? 4) 2 
iw zs) a 


< 
" 
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cg 7 127 fps 
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Therefore from (2): 


2v = 42.5 fps 


Therefore if no downwash effect is considered and accepting the 
assumptions that went into the momentum theory (i.e., constant v across 
disc, etc. ), the difference between the velocity in the wake and the free 
stream should be 42.5 fps. Assuming this the case, the engine was run 
statically and power adjusted to produce maximum wake velocities of 
40, 45, and 50 fps. These surveys were made with a velometer, which 
is a mass flow device. At these power settings wake surveys were made 
at the tail with the velometer and average tail velocities plotted versus 
wake velocity on Figure 20. 

Entering Figure 20 with wake velocity of 42.5 fps an average tail 
velocity of 31.2 fps was obtained. This velocity is then the prop wake 
contribution to tail velocity at the flight condition. From this tail effi- 


ciency was calculated as follows: 


st _ (84,5 + 31.2) 
t 4 (84.5)° 


. 88 


aS 
li 
= 


This figure is most idealized in that many assumptions were 
necessary to obtain it; however, it does justify values of a greater than 
1.0. If the downwash effect is considered an my = 1.3 is not considered 
to be out of order. Consequently, this value was used in the analytical 
development. 

It was also necessary to estimate the change in Cm; due to fuse- 


lage bending. The aircraft was suspended and a 25 lb weight placed on 
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the tail. At internal pressures of 5 psig and 7 psig the deflection of the 
tail relative to the wing root chord was found to be 1.8 degrees for 


5 psig and 0.9 degrees for 7 psig internal pressure. Therefore: 


Aa 
i 0.9 1 
—— = => = => deg/lb.f 
bet 7.0 25 Ceo 
Aa ) 
t 1 
— = —— deg/lb.f 
Net 5.0 13.9 


Then for any change in force on tail AF: 


AF 
mie. 2s 

AF 

oe a, 
Ger). pees *e 
or: 

AF 
Aa, = ae 


due to bending 


J+ 
pees ®, and q(8.46 psf at S. L.)a change in or of 1 degree will result 
in a change in force on the rigid tail of 8.4 Ib. f. 


Therefore: 


AB 76.4 Aq. 
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Now, the change in angle of attack of the tail due to the flexible fuselage 


may be expressed as: 


Aa 
t 
or: 
Aq 
t 
and; 
Aq 
t 
Solving for Aa ie 
Aa 
t 


Since C,.. 
at 
found previously: 


QO 
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= [ayer - Aa 


t 
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+ her = Aa 


rigid 
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mere: 


rigid 


varies directly with aor and substituting for values of K 
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From analytical development C,,,. = 1.397 
it ee. 
Therefore: zagad 
ot ) = 1.075 
t 7.0 
Ge. ) \e=oReni 
st 50 


These values will be compared with other methods of C,,. deter- 
i 


t 
mination in the results section of the report. 


Flight Calibrations 


The in-flight calibrations consisted of airspeed and normal ac- 
celeration calibrations. The cockpit airspeed indicator calibration was 
obtained by flying the aircraft at approximately 20 feet above the runway 
at indicated airspeeds of 45, 50, 55, 60, and 65 knots. These runs were 
made on reciprocal runway headings in order to compensate for wind. 
Knowing the runway length the average time for each run enabled a true 
airspeed calculation. This true airspeed was then corrected for tempera- 
ture and pressure altitude to obtain calibrated or corrected indicated air-~- 
speed. The results of the calibration are shown in Figure 21. Calibra- 
tions were run on January 3rd, January 7th, and April 25th with agree- 
ment on all three dates. The recorder airspeed indicator was calibrated 
against the cockpit indicator. However, as mentioned previously, the 
static source for the recorder airspeed indicator was a venturi mounted 
on the angle of attack boom mount. The aircraft pilot static tube was 
used as the dynamic pressure source for each instrument. 

The recorder accelerometer calibration was obtained by readings 


made of two cockpit accelerometers in steady turns. The aircraft was 
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equipped with a standard cockpit accelerometer and a glass tube accele- 
rometer was also added as a check. The two instruments were in agree - 
ment and the calibration of the recorder accelerometer did not change 


throughout the flight test. 


Flight Test of Stick Fixed Static Stability 


The static stability tests were conducted at a power setting to give 
54 knot IAS in level flight. (This was the condition used in the dynamic 
analysis.) This power setting was helid constant for all flight tests. 

An accurate record of fuel pitty was kept in order to know 
weight and center of gravity location for calculation of lift coefficients 
and correcting normal acceleration to a standard gross weight. Fuel 
was measured at the beginning and end of each flight. A consumption 


rate of 4 gallons (24 lbs) per hour was obtained and this rate was used 


to calculate fuel on board for each recorded event. 


Determination of Stick Fixed Neutral Points 


Neutral points were obtained for internai pressures of 7.0 and 
5.0 psig. The flight procedure was to set the power at 54 knots IAS in 
level flight. Elevator deflection and airspeed were then recorded at 
indicated airspeeds of 45, 50, 55, and 60 knots after a steady state 
flight condition had been reached. These runs were repeated by each 
pilot to obtain data at three different center of gravity locations. 

For the 7 psig runs an additional point was taken at 65 knots IAS, 
however, at 5 psig an undamped high frequency oscillation was en- 
countered at 63 knots IAS. An explanation of this motion is presented 
in the Discussion of Results. Figure 22 shows a typical trace for ele- 


vator versus airspeed and also a trace of the oscillation. 
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Elevator deflection was then plotted versus lift coefficient for 
each center of gravity location. These plots were linear and the slopes, 
oe Cis were plotted versus center of gravity location to obtain neutral 
point. Figures 23 through 29 show the plots of 6. versus Cr and 
Figures 30, 3+; and 32 show plots of é/ Cy versus xe for 5 psig, 7 psig, 


and rigid calculations. 


Determination of Stick Fixed Maneuver Points 


Maneuver points were also obtained for internal pressures of 
5.0 psig and 7.0 psig. For these test power was set at 54 knots JAS in 
level flight. The aircraft was then slowed to 45 knots IAS and pushed 
over to obtain 57 knots IAS. A pull up was then initiated with a constant 
elevator position. During the pull up airspeed, normal acceleration, and 
elevator position were recorded. The tests were repeated for each cen- 
ter of gravity location. At the point on the trace where the airspeed was 
54 knots IAS, acceleration and elevator positions were read and plotted 
on Figures 23 through 29. Figure 22 shows a typical trace of 6 V, and 
n. The normal acceleration was corrected toa standard gross weight of 


S/S 
420 lbs in order to match the flight condition used for the dynamic calcu- 


lations (i.e., ea otha gross weight, center of gravity at 29.00% MAC, 
and 54.0 knots IAS). The plots of be versus n were linear up to approxi- 
mately 2g's. At: load factors;sabove this the fuselage bending became non~ 
linear. However, since the dynamic motion was later found to be in the 
linear region, these linear slopes were used for plotting dé / dn versus 


ys to obtain maneuver points. Figures-30; 31, and 32 show plots of 


dé _/dn versus Xe, for 5 psig, 7 psig, and rigid calculations. 
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Determination of Longitudinal Stick Fixed Dynamic Stability 


Dynamic tests were conducted at 5 psig and 7 psig internal pres- 
sures. For the dynamic test the aircraft was loaded so as to obtain a gross 
weight of 520 lbs and a center of gravity location of 29.00% MAC in order to 
match the condition of the analytical development. 

The aircraft was excited with a pulse elevator input. It was found 
that some method of returning the pulse to the original amplitude and hold- 
ing it there was necessary. In order to achieve this a chord was attached 
to the pilot's lap belt and to a ring that could be slipped over the control 
stick. The length of the chord was adjusted to maintain 54 knots IAS in 
level flight. After each pull up the control stick was held against the 
chord to obtain a constant elevator deflection. This sort of 'Rube Goldberg" 
arrangement is typical of what must be resorted to when an airframe will 
not accept screws or nuts and bolts graciously. 

During the dynamic test elevator deflection, airspeed, angle of 
attack, pitch rate, and acceleration were recorded. Figure 33 is a typical 
recorder trace showing the dynamic response. 

In addition, during some of the tests the gun camera recorded 
fuselage deflection during the phugoid motion. A flutter of the elevator 
was used to correlate events on the camera and the recorder trace. From 
the recorder trace phugoid damping ratio (¢) was computed by Figure 18. 
Damping ratio and frequency of phugoid oscillation (@) were then plotted 
on Figure 17 which is a root locus plot obtained from the analytical 
development. 

The short period oscillation was not investigated. The analytical 
development indicated at short period time constant of the order of 0.2 
seconds. The angle of attack vane would not respond at this rate and from 


matching elevator input and angle of attack response it is impossible to 
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differentiate between aircraft short period and angle of attack vane response. 
In any event, the highly damped short period mode was not considered in- 
teresting enough to warrent a more intense investigation. 

Movie shots of fuselage bending were projected on a grid screen. 
From these observations it was possible to determine fuselage bending, in 
inches, at the tail and where in the motion the bending occurred. From the 
preflight test calibrations it was determined that the maximum bending for 
the 5 psig runs was of the order of 1.5 degrees, and 1.0 degrees at 7.0 
psig. This order of bending agrees with the magnitudes found in the pre- 
flight test determination of oon by statically bending the fuselage. Fig- 


ure 34 shows three typical frames from the movie shots. 


Determination of Cmi, and a 


In order to further verify the analytical development it was neces- 


sary to find the aircraft pitching moment due to tail incidence change oe 
and elevator effectiveness (7). 

To calculate Cm," flight test data from Figures 23 through 29 
was used to obtain values of on for center of gravity locations at a constant 
Cy of 0.5 (this was the value used for the rigid analysis). In this manner 
a slope me cee was obtained from Figure 35. 

From previous flight tests it was determined that the aircraft tail 
incidence was too high in that up elevator was necessary for almost ali 
flight conditions (note Figure 2). Therefore, after all other flight tests 
were completed the tail incidence was lowered 5 degrees. The aircraft 
was then flown at Cr of 0.5 at three center of gravity locations and elevator 
position recorded. As ee Ab. should not change with tail incidence a 
line was faired through these points with the a Ab. slope found pre- 


viously. 








oi 


From Reference 13: 
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and: Cms = ue Cmi, 
Therefore: 
Cy AX 
aes. = ee 
i, T 46 
e e 
and: 
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From Figure 35 D2 was found to be 0.55. Using this value of oe the slopes 


eee from Figure 35, and equation 1: 


Co = 61,10 
1 
ts 0 


Ci, = 71.23 
It 
7.0 


These values will be discussed in the results section of the report. 
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Other Aspects of Flight Test 


Some of the other unique features of the aircraft effected flight 
test procedures. The semi-open cockpit induced a chill factor that hada 
decided effect on pilot response during the late winter months. Also, the 
open cockpit proved a problem for voice transmission. This was solved 
by wearing an oxygen mask with valves removed. This procedure cut the 
background noise to a very low level. 

Due to the aft mounted propeller any loose gear in the cockpit 
could easily reach the propeller arc and be responsible for a fabric punc- 
ture. Thus, everything in the cockpit was of necessity tied down. 

The low aerodynamic mass of the aircraft made it very suscep- 
tible to light gust, and to obtain usable data extremely smooth air was 
required. Figure shows a noise trace ona day that would normally 
be considered quite smooth. 

Ground handling of the aircraft was reminiscent of sea plane 
operation and was a three man job. The ability of the aircraft to use 
sod runway greatly simplified the flight test, however. 

In all respects the three pilots concerned were quite pleased to 
be able to obtain a small taste of open cockpit and the scream of wires 


that was so typical of early aviation. 


ad 
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DISCUSSION OF RESULTS 


Definition of Problem 


The original purpose of this investigation was to determine the 
applicability of standard analytical techniques on predicting static and 
dynamic stability of an inflatable airplane. It was assumed that the 
major factor influencing any deviation from rigid analysis for this air- 
craft would be fuselage bending. A secondary effect would be wing tor- 
sion and bending affecting lift curve slope (ey) with changes with 
dynamic pressure (q). The change in Ci, with q was obtained from 
Reference 4 which is a report on full scale wind tunnel test of the air- 
craft. As this data was available and accounts for flexibility effects on 
the wing and tail, it was decided to concentrate on the effects of fuselage 
bending. 

The result of research was that there are many different opinions 
as to the effect of fuselage bending on such things as movement of stick 
fixed neutral and maneuver points (No and Nn’ It was then determined 
that a further investigation of the effects of fuselage bending on longi- 
tudinal stability might be profitable. This aircraft would seem to be an 
ideal vehicle for this type of investigation in that fuselage stiffness could 
be varied by changing internal pressure. 

All references agreed that fuselage bending was an overall de- 
stabilizing effect, and were in general agreement that Nn would move 
forward, except for Reference 9. Reference 9 stated that Nn would 
move aft due to a change in aircraft center of lift. The movement of 
Ny was not in agreement, however, in that there was difference of 
Opinion in how far ls should move forward relative to Nin 

It was anticipated that the change in Non and No between internal 


pressures of 5 psig and 7 psig would answer some of these questions. 
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It was considered that a valid measure of fuselage bending would 
be changes in C,,,, , and it was expected that Cm; would be reduced by 
i 


t 
fuselage flexibility. 


Predicted Results 


The rigid analysis was made using classical techniques of 
References 2 and 7. An effort was made to determine tail efficiency (n,) 
by a wake survey mentioned in the preflight test. This yielded a value 
of ne of 1.3 that resulted ina Cray, for the rigid aircraft of -1.397. A 
full scale wind tunnel value of Sis was chosen at the test q of 8. 46 psf. 

The results of this analysis were a predicted N, of 47.5 percent 
MAC, ye of 59.9 percent MAC, phugoid period of 15.7 seconds, anda 
damping ratio of 0.245. 

For the semi-rigid analysis Cray and Ci were factored out of 
the equations of motion and a root locus plotted showing effects of 
changing Cray: and “7 (Figures 16 and17). The excursion of Cr was 
expected to be from 3.0 to 5.0 and values of 3.0, 4.0, and 5.0 per radian 
were used for the root locus plot. 

A static bending test of the fuselage was made at internal pres- 
sures of 5.0 psig and 7.0 psig. (This procedure may be found in the 
flight test procedures section of the report.) The results of this test 


yielded predicted values of Cmi, of -1.08 at 7 psig and -0. 87 at 5.0 psig. 


Flight Test Results 


Although the flight test determination of Cm; and Ge were done 
at the end of the program a flight test value of T" tenecessary to corre- 
late other flight test results with predictions. From Figure 35, on was 
found to be 0.55. (The method used may be found in the flight test sec- 
tion of the report.) This value of we compared favorably with the pre- 


dicted value of 0.6 and is considered valid. 
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From the flight test for stick fixed static stability, N Ss 


O56 4 
found to be 40.0% MAC and Nor 9 42.0% MAC. Due to the scatter in 
data this result can only be considered as a qualitative indication in the 
movement of NG with fuselage bending. A possible explanation for the 
scatter in the neutral point data may be due to a change in position of 
the elevator potentiometer string tie down point. As velocity is 
changed, the leading edge of the horizontal stabilizer may deform and 
twist, consequently displacing the tie down point. Since the neutral 
point tests were made over a speed range of 45 to 65 knots, the tie 
down point would not remain in the same relative position, thereby 
making the elevator readings differ from the exact value in an unpre- 
dictable manner. It can be concluded that No moves forward with 
fuselage bending; however, the use of this data to conclude the rela~ 
tive movement of No and ie with flexibility is questionable. 


Although the slope of the od C, versus ae curve is not con- 


L 
sidered valid in finding Cm, , for a non-rigid aircraft, values of Sa, 
t 


were calculated from Figure 30. Cm; a found to be -1, 33 and 
hE 
Cm; , -1.54, using the flight test value of ae Using the data of 
a7. 
the 6 [Cy curves, Figures 23 through 29, as plotted on Figure 35, 
e 


was found to be -1.0 and Cm, 4 -1.23. This procedure 
tie 
is described in the flight test section of the report and is considered 


Ge 
a 6 


a much more valid method of determining Sioa, for a non-rigid air- 
craft than using the slope of the af C, versus ae curve. These 
values are in excellent agreement with those predicted by the static 
bending test. 

The flight test results for A were found to have much less 
scatter and are considered to be of significance. Nm5_ 9 was found 


to be 46.0% MAC and N 49.2% MAC. These results show a 


m7 4G 
definite destabilizing trend with bending. The slopes of the 6 /n 
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versus eee curves of Figure 31 yielded values of Ge of -1.07 and 
i 
Rieg 

of ~1.40. These values are in good agreement with magnitudes 


uaa 
predicted by the static bending test. 
Phugoid periods and damping ratios obtained from rigid analysis 
and flight test data were plotted on the root locus diagram, Figure 17. 
The results of these tests showed a definite trend towards instability 


with increasing fuselage flexibility. was found to be -0. 8 and 


Cmit 5.0 
Cmi, 2% -1.25 by the mean position of the plots of flight test periods 
and damping ratios at the two inflation pressures. 

The short period mode of dynamic motion could not be investi- 
gated with the instrumentation installed. It's high damping ratio of 0. 89 
and low time constant of 0.2 seconds precluded its usefulness for these 
tests. The reason for the high damping ratio is the very low aerodynamic 
mass of this aircraft. 

Table III is a comparison of Sm values found from static bending 
tests and four flight test methods. This comparison shows close agree- 
ment by all methods and a definite destabilizing trend with fuselage bend- 
ing. 

The results of flight test thus verified the values of Cmi, and o. 
used in the rigid analysis. It can then be stated with some certainty 
that the rigid analysis was valid. This being the case, Ny moved for- 
ward from the rigid value to the 5 psig value 7.5% MAC while No moved 
forward 13.9% MAC. While the movement of Ne between 5 psig and 
7 psig is not well defined, its movement from a rigid value to a flexible 
value is considered to be of such validity as to warrant mention. 

It may therefore be stated that for this aircraft es and Ne are 
both destabilized by fuselage flexibility and it is further indicated that 
ee moves further forward than No’ 

The high frequency aircraft oscillation mentioned in the flight 


test procedures was determined to be caused by elevator flutter. This 
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flutter was found to become excited at 63 knots IAS with 5 psig internal 
airframe pressure. 

A movie supplement to Reference 4 relates a similar motion 
occurring at 7 psig pressure at a dynamic pressure of 20 psf in the 
full scale wind tunnel. The motion was felt to be excited by elevator 
flutter in that fixing the elevator to the horizontal stabilizer stopped 
the motion. 

Figure 22 shows a recording of the motion during flight test at 
5 psig internal pressure. From the recording it was found that the ele- 
vator oscillation was from amplitudes of + 2 degrees to + 13 degrees at 
a frequency of 5 cps. It is felt that this is the same motion described 
in the wind tunnel test. 

In the wind tunnel analysis the cause of the flutter was attributed 
to control looseness. This assumption is felt to be valid in that the 
elevator can be moved approximately 5 degrees with the stick fixed at 
7 psig internal pressure. Ata lower pressure the elevator is, of 
course, more loose. It is felt that, this being the case, the elevator 
flutter could occur at a dynamic pressure of 12 psf (63 knots IAS), and 
that the motion observed in flight test is the same motion previously 


seen in the wind tunnel. 
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CONCLUSIONS 


It is concluded that rigid analysis techniques will predict longi- 
tudinal stick fixed dynamic stability of the YAO-3G inflatable airplane 
to a good degree of accuracy. Static stability predictions are not as 
accurate but do indicate typical values. The aircraft when operated at 
recommended internal pressures is of sufficient rigidity as not to 
deviate far from rigid analysis. However, if rigid values of Cmi, are 
corrected by factors obtained from static bending tests, accurate pre- 
dictions may be obtained for both dynamic and static stick fixed sta- 
bility. 

The effect of fuselage flexibility of this aircraft was destabiliz - 
ing to both dynamic and static stability. The phugoid mode was ren- 
dered less stable by increased flexibility. Both stick fixed neutral 
and maneuver points moved forward; however, it was indicated that 


Nun is destabilized toa larger extent than N . 
° 
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TABLE I 


Aircraft Specifications 


General: 
Span PAA TS ai 
Length 22, 0it 
Height Tee wee 
Wing area Zane 
Wing loading Ae 
Thrust loading 1.69 #/ £t* 
Moment of Inertia (ly) 277.2 slug ft? 
Weights: 
Weight empty 315 lbs 
Useful load - 1 pilot 194 lbs 
Fuel 10 gal 61 lbs 
Gross take-off 570 lbs 
Engine: 
Available sea level static thrust 220 lbs 
at combat rating at 3750 RPM 
Propeller diameter 47 in 
Height of thrust line above CG 3.1446 


Minimum fuel consumption at 


sea level cruise 20.2 #/hr 
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TABLE I (Continued) 


Airfoil section NACA 0015 
Area (S) N24 Sete 
Span (b) 23a. ft 
Aspect ratio (AR) 4.44 
Taper ratio (A) 1 
Chord (c) 53) it 
Ailerons: 
Chord ci /c 20.4 % 
Area Sts 6.41 ft® 
Horizontal Tail: 
Area (S ) ZIOb tt 
Span (b,) 6.66 ft 
Chord (c.) S, a SaaliGuet’ 
Tail volume NE ae ~399 
Elevators: 
Chord (c.) 1/33 4t 
Area 8.86 ft? 
Elevator travel 30 up to 25 down 


Vertical Tail: 


Area (S_) lomo fee 
“Vv 
Span (b) SO Ceare 
Chord (c ) 4.416 ft 
Vv 
Rudder: 
Mean chord 1.66 ft 


Area 5.0 ft? 
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TABLE II 


Analog Computer Potentiometer Values 


Quantity Value Pot Setting 
1 -1 


mab 








soe 
. 2418 


2338 | 


. 2502 


C= 
™it/ 2 
. 681 
. 3625 


= 40.3 


K 


Note: 


Voltage-quantity ratios 
da, A6,.45.: Iv. = .01 rad Au, Av: lv. = 1 ft/sec 


Memo 1 y. = Ol aes sec Ah: lv. = 1 ft 
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TABLE Il 


Cm; Values 


Rigid Analysis using 


Classical Techniques 


Static Bending Test 


(with factors applied 
to rigid values) 


Static Flight Tests 


a) gradient of 


Dynamic Flight Tests 
Figure 17 














Three View Drawing 


Figure |. 
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Figure 33. SFIM Recorder Tapes 








Level Flight Condition 
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Up Load 
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Down Load 





Fig. 34. Selected Prints From Movie Camera Showing Fuselage Bending. 
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